Transport in gapped bilayer graphene: the role of potential fluctuations 
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We employ a dual-gated geometry to control the band gap A in bilayer graphene and study the 
temperature dependence of the resistance at the charge neutrality point, Rnp(T), from 220 to 1.5 
K. Above 5 K, Rnp(T) is dominated by two thermally activated processes in different temperature 
regimes and exhibits exp(Ts/T) 1//3 below 5 K. We develop a simple model to account for the ex- 
perimental observations, which highlights the crucial role of localized states produced by potential 
fluctuations. The high temperature conduction is attributed to thermal activation to the mobility 
edge. The activation energy approaches A/2 at large band gap. At intermediate and low tem- 
peratures, the dominant conduction mechanisms are nearest neighbor hopping and variable-range 
hopping through localized states. Our systematic study provides a coherent understanding of trans- 
port in gapped bilayer graphene. 



PACS numbers: 72.80.Vp, 73.22.Pr, 72.20.Ee 

Bilayer graphene is a unique two-dimensional (2D) ma- 
terial with a tunable bandgap. A perpendicular elec- 
tric field breaks the inversion symmetry between the 
two graphene layers and results in a field-dependent 
bandgap 1 3 . Its experimental signatures have been ob- 
served by infrared spectroscopy^' and angle-resolved 
photoemission 8 , but remain incomplete and perplexing 
in transport!^ 2 -!. Near room temperature, Xia et al. ob- 
serves thermally activated conduction and attributes it 
to Schottky barriers at the electrode-gapped bilayer in- 
terface 11 . In the mK regime, Oostinga et al. reports 
variable-range hopping 10 . To date, systematic investiga- 
tions combining high and low temperatures are lacking 
and a coherent understanding of conduction in gapped 
bilayer has yet to emerge. 

In this work, we control the band gap in bilayer 
graphene using top and bottom gates, and measure the 
temperature-dependent resistance at the charge neutral- 
ity point (CNP) Rnp(T) as a function of the gap, in 
the temperature range of 1.5 K < T < 220 K. We de- 
velop a model to explain the data, which highlights the 
essential role of localized states produced by potential 
fluctuations. Our data point to three conduction mech- 
anisms: thermal activation to the mobility edge at high 
temperatures, nearest neighbor hopping at intermediate 
temperatures and variable-range hopping at low temper- 
atures. 

We fabricate Si02/Hf02 dual-gated bilayer graphene 
field effect transistors using procedures previously de- 
scribed in Ref. [T3l 30 nm Hf02 is deposited on single 
or bilayer graphene by atomic layer deposition and used 
as the topgate dielectrics. We have achieved high mobil- 
ity \i of 9,000 - 16,000 cm 2 /Vs on single layer graphene^. 
Here, on dual-gated bilayer, \i ranges from 1,500 to 6,000 
cm 2 /Vs, which is generally lower than /i up to 12,000 
cm 2 /Vs observed in our pristine bilayer samples. Ra- 
man spectra on dual-gated bilayer devices show no visible 
D band, indicating minimal defect creation (Fig. S4 in 
Ref. |23j) . The gating efficiency of the topgate is approxi- 
mately 2.8xl0 12 /cm 2 per volt from Hall measurements, 
which is ~ 40 times of the efficiency of the 290 nm Si02 




FIG. 1: (Color online) R(Vt s ) at fixed Vb g of sample Bl at 
10 K. From left to right, V hg = +30 V to -60 V at 10 V steps. 
Inset: a SEM picture of a device. The mobility \i is approx- 
imately 1,500 cm 2 /Vs for sample Al, A2 and 3,000 cm 2 /Vs 
for sample Bl, B2. 

backgate. The gating range of the topgate is greater than 
1.4xl0 13 /cm 2 . 

A scanning electron microscope (SEM) micrograph of 
a dual-gated device is shown as an inset to Fig. [I] Two 
types of topgates are used. In samples Al and A2, the 
area between the two voltage probes is completely cov- 
ered by the topgate electrode (top of inset), so that it 
does not include the interfacial resistances between the 
gapped and ungapped area. In samples Bl and B2, the 
coverage is partial (bottom of inset). All four measure- 
ments yield similar behavior, indicating that the resis- 
tance of gapped bilayer is dominated by the bulk. 

To obtain the T-dependence of the charge neutrality 
point resistance R^p(T), we fix the backgate voltage Vb g 
and sweep the topgate V tg . The sweep is repeated at 
different Vb g 's and temperatures. Figure [I] shows i?(Vt g ) 
of sample Bl at 10 K. From left to right, Vb g changes 
from +30 V to -60 V at 10 V steps. The maximum re- 
sistance of each curve corresponds to the CNP for that 
specific pair of (Vbg? Vt g ) settings. The minimum of all 
i?NP is found at roughly Vb g o = V and V^ g o = 3.1 V 
in Fig. [I] which corresponds to the condition where the 
average band gap A = G 4 . The offsets Vt g o and Vb g o 
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FIG. 2: (Color online) (a) R^p{T) of sample Bl in a semi-log 
vs. l/T plot in the temperature range of 1.5 K < T < 220 K. 
From top to bottom, V^ g = -60 V to -10 V at 10 V steps. The 
triangles mark 50 K and 5 K respectively, (b) The expanded 
region of 70 K< T < 220 K in (a). The solid and dashed lines 
in (a) and (b) are guide to the eye. 

comes from unintentional chemical doping. At the CNP, 
A was shown to increase with the electric displacement 
field D = D hg = —D tg = e hg V^ g /d hg , where e bg and d hg 
are the dielectric constant and thickness of the backgate 
oxide respectively, and V£ = Vbg — VbgO is the effec- 
tive backgate voltage^. Using esi02 = 3.9 and 6fei02 = 
290 nm, we obtain D = 1.34 V/nm for V£ g = 100 V. 
In samples Bl and B2, the resistance of the dual-gated 
area 2 is calculated by subtracting the resistance of ar- 
eas 1 and 3 (Fig. [I] inset), which is obtained by sweeping 
Vb g while grounding V tg at each temperature. Standard 
lock-in techniques are used with excitation current 0.2 
nA to 50 nA. Low currents are carefully chosen for high- 
resistance measurements to avoid current heating. 

Rnp(T) of sample Bl is shown in a semi-log vs. l/T 
plot in Fig. [2|a). From top to bottom, V£ changes from 
-60 V to -10 V at 10 V steps. R^p(T) shows an insulat- 
ing T-dependence in the whole temperature range. The 
high-T region is expanded in Fig. [2^b), where logi?NP vs 
l/T shows a linear trend (solid lines) at T > 150 K. This 
trend suggests thermally actived conduction. From 150 
K to 50 K, the slopes continue to decrease, and become 
approximately constant again between 50 K and 5 K al- 
beit with much smaller values (dashed lines in Fig.[2|a)). 
This suggests another thermally activated process with a 
much smaller activation energy. Below 5 K, R^p(T) re- 
mains nearly T-independent for small | and increases 
slowly with decreasing T for larger \V£ \. Similar behav- 
ior is observed on all four samples. 

These observations are remi nisce nt of impurity conduc- 
tion in doped semiconductor d 14 l 15 l, where localized states 
produced by impurities reside in the band tail and are 
separated from delocalized states by a mobility edge E c . 
In such systems, thermal activation to the mobility edge 
dominates the high-temperature conduction. As T de- 
creases, electrons hop through localized states, overcom- 
ing the energy difference between adjacent neighbors, i. 
e., nearest neighbor hopping. At even lower T, elec- 
trons hop through localized states far from each other 
via the assistance of phonons of energy /crT, giving rise 
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FIG. 3: (Color online) R^p(T) of sample Al in the tempera- 
ture range of 1.5 K < T < 220 K. From top to bottom, V£ s 
= -90 V to -30 V at 10 V steps. The solid lines are fittings to 
Eq. ^ The 3 terms in Eq. ^are plotted as dash-dotted, dot- 
ted, and dashed lines respectively for the top curve to show 
the contribution of each term at different temperatures. 

to variable-range hoppingp^. 

Localized states due to disorder also exist in bilayer 
graphene^H Scanning tunneling microscope (STM) mea- 
surements show the potential of the CNP in bilayer can 
fluctuate between ± 40 meV 17 . Sim ilar ob servations have 
been made on single-layers as wel l 18 * 19 *. The resulting 
electron and hole puddles are typically a few tens of nm 
in lateral siz d 17 " 19 i In a gapped bilayer, these puddles are 
confined and separated from the mobility edge above the 
potential fluctuations by a binding energy. In a zero-gap 
bilayer, Klein tunneling may provide a p aralle l conduc- 
tion path to thermally activated processe d 20 * 21 !. 

The T-dependence of R^p(T) in Fig. [2] and the above 
reasoning led us to propose a fitting including the follow- 
ing three terms: 

^np(T)- 1 = R^ 1 exp[-E 1 /k Bj T] 
+R2 1 exp[-E 2 /k B T} 
+ R-iexp[-(T 3 /T) 1 /*] (1) 

where E\ , E 2 and T 3 represent the two activation and the 
hopping energy respectively and R 2 and Rs are the 
corresponding resistance coefficients. Eq. [I] produces an 
excellent description of i?NP (T) in the whole temperature 
range in all samples. Each parameter plays a different 
role in the fitting 2 ^. The fitting curves of sample Al are 
shown in Fig. [3] The combination of different energy 
scales accurately captures the shoulder seen around 50 
K (Fig. [3J, and the bend around 5 K (Fig. ^a)) in all 
samples. 

The extracted E\ of all four samples are plotted in 
Fig. |4ja) and show similar trends. E\ remains approxi- 
mately constant at 250 K for small |V^ g | (or equivalently 
small D), and starts to increase at larger |V^ g |, reaching 
- 450 K at VI = -90 V in samples Al and A2. Also 
plotted in Fig. 0a) is the calculated A/2 of Ref. 1 (solid 
black line). At small |V^ g |, E\ differs from the small 
A/2 predicted by theory but approaches theory with in- 
creasing | V^g I- I n particular, a finite E\ at V£ = V, 
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FIG. 4: (Color online) (a)(b)(c) E u E 2 , and T 3 vs. V b ' g for 
sample Al (solid black square), A2 (hollow red square), Bl 
(solid blue triangle) and B2 (hollow green triangle). The the- 
oretical A/2^is plotted as a solid line in (a) for comparison. 



where A = 0, suggests a different activation mechanism. 
We attribute E\ to the activation barrier to the mobility 
edge, i. e., E\ = E c — Ep. Since at A = 0, the Fermi 
level Ey lies approximately in the middle of the fluctu- 
ating disorder potential <£, we associate E\ = E c with 
the root mean square (RMS) amplitude of the potential 
fluctuation, i. e., $rms = E\ = 250 K or 21.5 meV. 
Independently, <£rms m the STM data of Ref. [17] was 
found to be 19 meV^, in good agreement with results 
obtained here. In this scenario, all carriers residing in the 
puddles will be excited above E c at very high tempera- 
tures. Estimating the carrier density in the puddles to 
be n = $ R Ms(2ra77rft 2 ) = 6x 10 n /cm 2 (ra* 0.033ra e ) 
and using /i = 1,500 - 3,000 cm 2 /Vs for all samples, we 
expect Ri = 1/ne/i to be 3500 - 7000 Q/\3. R\ extracted 
from the fittings ranges 3,000 - 5,000 fi/D, in excellent 
agreement with our estimates^. 

The above analysis should also hold for a small band 
gap A < $rms? which explains the nearly flat E\ at small 
|Vb g | in Fig.j4|a). As A further increases, we expect E\ to 
be gradually dominated by the band gap itself and even- 
tually approach A/2. Our results in Fig. [4^ a) support 
this trend. Small differences likely arise from the details 
of the disorder potential in individual samples. The 
above model also explains the insulating T-dependence 
widely observed near the CNP in backgated bilayers^HUl 

As T decreases, the above thermally activated conduc- 
tion becomes inefficient; our data in Fig. [2] show that 
other low-energy processes gradually take over from 150 
K to 50 K. Between 50 K > T > 5 K, the rise of R K 



(T) is largely given by the 2 nd exponental term in Eq 
(Fig. J3J) . The extracted activation energy E2 is shown in 
Fig. Hfb). E 2 is nearly V^-independent and averages ~ 
10 ICJor 25-50 times smaller than E\, in all samples. In 
doped semiconductors, E2 is associated with the hopping 
conduction between nearest neighboring impurity states. 
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FIG. 5: (Color online) (a) A simple model of the potential 
fluctuations in bilayer. The bound state of each cylindrical 
well is schematically shown. A wide well may have more than 
one bound state. (b)The bound state energy £i(L) for Vb = 
43 meV (solid magenta line) and Vb = 13 meV (dashed black 
line) . 

Percolation theory gives E2 = < > = < > is the 
average energy difference between neighboring localized 
states on the percolation path 14 . Similar to the localized 
states in doped semiconductors, in a gapped bilayer, elec- 
tron and hole states near the band edges are quantized 
and localized due to small puddle size 2 ^. These localized 
states may support a hopping conduction that is more 
effecitive than the thermal activation above the mobil- 
ity edge at k^T « E\. In this scenario, the hopping 
is expected to occur between nearest neighbors (NNH) 
when the temperature is not too low, and transitions 
to variable-range hopping (VRH) at the lowest tempera- 
tures^. These two hopping processes can account for the 
exp[— E2/k^T] and the exp[— (T3/T) 1 / 3 ] terms necessary 
to describe our data. 

The following analysis provides an order of magnitude 
estimate of E 2 = < > in our samples. We assume 
the distribution of the puddle size L in our samples is 
similar to that in Ref. 17, since both samples are on Si02 
substrates and have similar $rms- This distribution is 
estimated by analyzing the potential map of the CNP in 
Ref. [171 A histogram of L is generated and modeled by a 
Gaussian distribution as shown in Fig. S2(b) 23 . The best 
fit yields L = (13±8.4) nm, consistent with other STM 
reports^. In our simplified model, the disorder potential 

is approximated by a random network of cylindrical 
wells with a single depth Vb but varying diameter L as 
shown in Fig. [5|a), neglecting the shape variation and 
the distribution of Vb- The bound state energy e\ of each 
well is calculated using a constant effective mass m* = 
0.033raJ^. Fig.jSJb) plots s- x vs. L for two representative 
Vo's. The solid curve corresponds to Vb = 43 meV, which 
is 2<£rms in our samples. The dashed curve corresponds 
to Vb = 13 meV. We obtain the average < > using 
£i (L) and the distribution of Zp3 including puddles of all 
S i ze p] We find < £[ . > = 11Q K and 9 K for V = 43 

meV and 13 meV respectively. In real samples Vb obeys a 
Gaussian distribution 23 . Statistically 25% of all puddles 
have Vb < 13 meV and 68% of all puddles obey Vb < 43 
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meV. Thus these calculations provide reasonable bounds 
for < £ij > in our samples. The above analysis suggests 
that < £ij > is likely to be a few tens of Kelvin, which is 
in reasonable agreement with our extracted E 2 ~ 10 K 
and supports the NNH scenario. 

The above model points to the onset of VRH conduc- 
tion at yet lower temperature. Indeed, R^p(T) data in 
Fig. [2|a) shows a bend at 5 K. We have verified that 
this bend is not due to current heating. Although the 
temperature range is limited, the VRH term in Eq. [I] is 
necessary to capture this behavior. At large |T^ g |, the 
extracted T3 ranges 0.7 K - 1.3 K as shown in Fig. gc). 
These values are in very good agreement with T 3 = 0.5 
- 0.8 K obtained in Ref. [TUJ where the VRH conduction 
is observed in similar samples for a larger temperature 
range 5 K > T > 50 mK. At small \V^ g \, T 3 is so small 
that Rnp(T) remains essentially T-independent, proba- 
bly due to Klein tunneling. 

The corresponding R3 is ~8 kfi/D at small |Vj£ | and 
increases with increasing \V£ \ in all samples. R% reaches 
approximately 75 kO/D in samples Al and A2 at |V^ g | 
= 90 V. This trend is consistent with the increasing sup- 
pression of Klein tunneling as A increases. The opening 
of the gap may also alter the strength of the electron- 
acoustic phonon coupling, which is essential to the VRH 
process, and affect R%. As a final check, we note that 
VRH is expected to occur at T « T c , where T c is de- 
termined by Eq. |2j where the variable activation energy 
E a becomes much smaller than E^^. 

Ea{Tc) - a(k B T)-i lT = T «-— T' (2) 



Using T 3 = 0.9 K and E 2 = 10 K, we obtain T c = 33 
K. Thus, the assignment of VRH at T < 5 K is self- 
consistent. 

In conclusion, we measure the T-dependent resistance 
at the charge neutrality point in gapped bilayer and de- 
velop a simple model to account for the observed multi- 
ple energy scales. Localized states produced by potential 
fluctuations play a crucial role in our model. At high 
temperatures, thermal activation to the mobility edge 
dominates the conduction. The activation energy is de- 
termined by the potential fluctuations at small bangap 
and approaches that of a band insulator at large bandgap 
or in clean samples. At lower temperatures, a percola- 
tion network forms via carrier hopping through localized 
states, leading to a second activated process and eventu- 
ally variable range hopping. Our model provides impor- 
tant insights into electrical transport in gapped bilayer 
graphene, which may be useful in a range of electronic 
and optical applications. 
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The states are localized when the bandgap is sufficiently 
large so that the puddles are separated by tunnel barriers. 
When the gap is small, Klein tunneling should take place 



although calculations show that it may be suppressed by 
strong disorder^. Experimentally, the weak T-dependence 
of .Rnp (T) at low Vb g renders the distinction of the 2 nd 
and 3 rd term in Eq. (1) less meaningful. These subtleties 
should be kept in mind when looking at E2 extracted for 
small Vbg- 

Ref. [3 predicts a much larger m* near the CNP in gapped 
bilayer. A heavier mass means lower Si and more bound 
states in a well, which may lead to lower < £ij >. This 
effect has not been included in our simulations. 
Puddles on the percolation path may have a narrower en- 
ergy distribution and consequently lower < >. This 
effect has not been considered in our simulations. 



